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ABSTRACT 


The  finite-difference  method  is  used  for  the  nonlinear  analysis  of 
shells  of  revolution  consisting  of  elastic  shell  segments  of  various 
geometries  and  wall  constructions  joined  by  elastic  rings.  The  analysis 
and  associated  digital  computer  program  were  developed  in  response  to 
the  need  for  a  better  design  tool  for  practical  shell  structures. 
Numerical  results  are  presented  for  displacement  and  stress  distributions 
in  various  pressure  vessels.  Particular  emphasis  is  given  to  systems  in 
which  nonlinear  effects  are  important  and  may  influence  the  design. 

Values  calculated  with  linear  theory  are  compared  with  those  from  non¬ 
linear  theory. 
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Introduction 


The  work  described  herein  represents  part  of  an  effort  to  develop  a  user-oriented 
computer  program  for  the  analysis  of  the  stress,  stability,  and  vibration  of  elastic, 
ring-stiffened,  segmented  shells  of  revolution  with  various  wall  constructions  and 
submitted  to  various  axisymmetric  load  systems. 

Such  a  structure  is  shown  in  Fig.  1,  for  example.  This  shell  consists  of  four 
segments,  labeled  (T)  ,  (2)  ,  (5)  ,  and  (?)  .  Segment  (l)  is  a  meridionally 
stiffened  coverplate  or  hatch;  (2)  is  a  spherical  segment  reinforced  by  T-shaped 
rings;  (5)  is  a  cylinder  with  ring  stiffeners;  and  @  is  an  ellipsoidal  closure  of 
variable  thickness.  The  thickness  of  all  of  the  segments  varies  near  the  junctures. 
While  Fig.  1  does  not  represent  an  actual  design,  it  does  illustrate  some  of  the  features 
incorporated  into  the  design  of  practical  shell  structures.  The  shell  structure  may  be 
loaded  by  internal  or  external  pressure,  and  line  loads  and  moments  may  be  applied 
to  any  or  all  of  the  ring  reinforcements.  The  designer  may  be  interested  in  the 
behavior  of  such  a  composite  shell  structure  with  respect  to  stress,  stability,  and 
vibration. 

A  user-oriented  computer  program  called  BOSOR  2  has  been  written  for  the 
general  analysis  of  shells  of  the  type  shown  in  Fig.  1.  This  program  is  a  "next 
generation"  successor  to  the  code  BOSOR  described  in  Refs.  1  and  2.  The  program 
BOSOR  2  consists  of  two  parts:  a  part  in  which  the  axisymmetric  stresses  and  dis¬ 
placements  are  calculated  from  nonlinear  theory,  and  a  part  in  which  nonsym  metric 
loads  and  vibration  frequencies  are  calculated  from  a  linear  eigenvalue  analysis.  The 
linear  stability  and  vibration  analyses  with  numerical  results  are  given  in  Ref.  3. 
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In  this  paper  the  nonlinear  analysis  for  axisymmetric  stresses  and  displacements  is 
described,  and  numerical  results  are  given  for  examples  in  which  nonlinear  effects 
are  significant. 

The  method  of  finite-differences  is  used  to  calculate  axisymmetric  stress 
states  in  segmented,  ring -stiffened,  elastic  shells  of  revolution  with  various  types 
of  wall  construction.  Nonlinear  equations  similar  to  those  developed  by  Rcissner 
(Ref.  4)  are  used  in  the  analysis.  The  work  is  an  extension  of  that  presented  in 
Ref.  5.  The  applicability  of  the  analysis  of  Ref.  5  is  extended  to  include: 

1.  Analysis  of  segmented  (composite)  shells,  such  as  cylinder-cone 
combinations  or  joined  shells  with  dissimilar  wall  characteristics. 

2.  Analysis  of  shells  with  discrete  rings  at  a  number  of  stations  along  the 
meridian  and  at  the  boundaries. 

3.  Analysis  of  shells,  the  wall  properties  of  which  vary  in  the  meridional 
direction. 

4.  Increased  generality  of  the  type  of  axisymmetric  loading  applied  to  the 
shell. 

In  addition,  much  effort  has  been  devoted  to  making  the  computer  program 
user-oriented.  It  is  felt  that  the  main  contribution  of  this  paper  to  the  art  of 
shell  analysis  is  the  development  of  a  design  tool  for  practical  shell  structures 
with  ring  supports,  in  which  nonlinear  effects  may  be  important.  In  this  respect 
the  paper  represents  extensions  to  the  work  of  Kalnins  (Ref.  6),  Radkowski, 
et  al  (Ref.  7),  Mason,  et  al  (Ref.  8),  Wilson  and  Spier  (Ref.  9),  and  Scpetoski, 
et  al  (Ref.  10). 
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One  of  the  more  important  extensions  of  the  analysis  of  Ref.  5  is  the  addition 
of  a  capability  to  treat  segmented  shells.  This  capability  can  also  be  used  to  advantage 
for  simple  shells.  For  instance  a  rather  long  cylinder  submitted  to  pressure  loading 
can  be  divided  artificially  into  three  segments :  two  edge  segments  in  which  stresses 
and  displacements  vary  rapidly  over  short  lengths,  and  a  central  segment  in  which 
deformations  are  uniform.  The  station  spacing  in  the  finite  difference  mesh  can  be 
small  in  the  two  edge  segments  and  large  in  the  central  segment.  It  has  been  found 
from  experience  that  it  is  better  to  divide  the  shell  into  segments,  and  thus  to  maintain 
uniform  station  spacing  within  each  segment  rather  than  to  vary  the  station  spacing 
within  any  segment.  Examples  are  given  in  the  section  on  Numerical  Results  in  which 
single  shells  are  treated  in  segments. 


Nomenclature 


AU«  Ai2’  etc* 

E 

e 

H 

•1 

H 

I.I.I 
x*  y*  xy 

i 

K 

M 

M 


ring  cross-section  area 

coefficients  of  constitutive  equations  (6) 

Young's  modulus 

ring  shear  center  eccentricity,  see  Fig.  3 
radial  load  applied  to  ring,  see  Fig.  3 
shell  radial  load,  see  Fig.  2 

moments  of  inertia  about  x,  y  axes;  product  of  inertia  (see  Fig.  3) 
mesh  points 

Gaussian  curvature  1/R.  R_  ;  also  total  number  of  mesh  points 
couple  applied  to  ring,  see  Fig.  3 
shell  moment/length,  see  Fig.  2 
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N 


P 

r 

R 

u 

V 

V 
P 
e 

K 

<P 

<t> 


stress  resultants,  see  Fig.  2 

pressure,  see  Fig.  2 

parallel  circle  radius,  see  Fig.  2 

Radius  of  curvature 

radial  deflection,  see  Fig.  2 

axial  load/length  in  shell,  see  Fig.  2 
axial  load/length  applied  to  ring,  see  Fig.  3 

meridional  rotation,  see  Fig.  3 

reference  surface  strain 

reference  surface  change  of  curvature 

angle  from  horizontal  to  tangent  to  deformed  meridian,  see  Fig. 
stress  function  ip  =  rH 


Subscripts  and  Superscripts 

.n 


(  > 
<>+ 
(  f 

<>s 

<>r 
(  )’ 

(>v 

<>H 

<>2 

a 


nth  iteration 

after  meridional  discontinuity,  see  Fig.  3 
before  meridional  discontinuity,  see  Fig.  3 
shear  center  of  ring;  also  "shell" 
pertaining  to  ring 

derivative  with  respect  to  arc  length 

pertaining  to  axial  (vertical) 

pertaining  to  radial  (horizontal) 

meridional 

circumferential 

undeformed 


LOCKHEED  PALO  ALTO  RESEARCH  LABORATORY 

lOCIHIIt  Mltllltt  *  *  P  A  C  t  COMPANY 

A  OIOUP  aivlllOM  O*  lOCPMItO  .ItCtAPt  CORPORATION 


Basic  Equations 


The  nonlinear  equations  of  equilibrium  of  a  shell  of  revolution  subjected  to 
axisymmetric  loading  are  presented  in  Ref.  4.  They  are: 


(rV)  +  rpv  *  0  -  N2  +  rpH  =•  0 


(rMj)1  -  cos  +  ^  8*n  *P  ~  008  <p  =  0 

where  the  sign  conventions  are  shown  in  Fig.  2.  Here  tp  -  rH  and  superscript  prime 
indicates  differentiation  with  respect  to  arc  length  s  .  The  strain-displacement  and 
curvature -displacement  relations  given  in  Ref.  4  are: 


cos  <p  /  i  \ 

€i  =  i1  +  -  1 


e2  =  u/r 


“l  =  <  - 


«2  =  (sin  -  sin  <p)/ r 


and  the  compatibility  equation  derived  by  elimination  of  u  from  Eqs.  (2a)  and  (2b)  is 


cos  <PQ(re2) '  -  cos  (pr'e^  =  r'(cos  <p  -  cos 


The  basic  equations  are  the  compatibility  Eq.  (3)  and  the  third  equilibrium  Eq.  (lc). 
The  basic  dependent  variables  are  <p  and  4>  .  The  strains  and  their  derivatives  arc 
assumed  small  compared  to  unity.  With  r1  =  cos  <pQ  Eqs.  (lc)  and  (3)  can  be  written 
in  the  form 


cos  <px  =  f 
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where  in  Eq.  (lc): 


in  Eq.  (3) 


X1 -  M1 

X1  =  e2 

x2  =  M2 

X2  =  £1 

(5) 

f  =  rV  cos  <p  -  ip  sin  <p 

f  =  cos  (p  -  COS 

The  constitutive  equations  for  orthotropic  shells  can  be  written  in  the  form 


V 

A11  A12  A13  A14 

'Nl' 

> 

> 

> 

> 

N„ 

2 

21  22  23  24 

2 

M1 

A31  A32  A33  A34 

K1 

M2. 

A41  A42  A43  A44. 

' 

K 

2 

Reissner  (Ref.  4)  formulated  the  equilibrium  and  compatibility  equations  for  the 


case  of  moderate  rotations  by  assuming  that 


cos  <p  =  cos  {q>Q  -  P)  «  (l  -  |/32)  r'  +  0r/R2 
sin  <p  =  sin  [<pQ  -  p)  «  (l  -  %P2)*/R2  ~  /3r' 


(7) 


The  dependent  variables  of  the  nonlinear  stress  analysis  are  now  p  and  V'  •  Since 

the  equilibrium  and  compatibility  equations  (4)  are  written  in  terms  of  c  ,  c  ,  M  , 

1  z  1 

and  M2  ,  which  through  Eqs.  (6)  are  expressed  in  terms  of  ,  N2  ,  ,  and  « 2  , 

one  must  obtain  the  latter  quantities  in  terms  of  p  and  ip  .  thus: 

rfcj  =  tfi  cos  tp  +  rV  sin  <p 

N„  =  ip'  +  rp„ 

2  H  (8) 
=  p' 

*2  =  (pr>  +  i  02r/R2)/r 
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Substitution  of  the  right -hand-sides  of  Eqs.  (8)  in  Eqs.  (6)  and  of  Eqs.  (6c,  d)  in  the 
left -hand -side  of  Eq.  (4)  and  use  of  Eqs.  (7)  for  cos  <p  and  sin  ip ,  leads  to: 


rA32*"  +  lr'(A31  -  A^  +  A32)  +  rA^]  *'  +  [r'A^  -  r'2A41/r  -  A31rK]  *  +  rA^i" 
*  lr'(A33  +  A34  -  A43  -  rVA,,)  +  rA;,]/3'  +  [r'A^  -  r,2A^/r  -  A^rK 


31'  33 


+  rV(AgjrK  -  rAjj  +  r  A4j/r  -  A4^r/R2)  +  A3Jrr  py  -  A42r  p^/R^lP 
+  (A31  ”  A42)(r/R2)  4>'P  +  (A31r,/R1  +  rA^/R2  -  2r'A41/R'2)^  +  (rAg/Rg) 
~  <A43  -  A34)(r/R2)^'  -  |  (3A44r’/R2  -  A^r'/Rj  -  A^r/R^p2 
+  rV  <A31rVRl  +  A'ir/R2  -  A41r'/R2)  +  A32  (r\)'  +  r2pH(A^  -  r'A^/r) 


"  A31r  pt/R2 


(9) 


The  expression  (9)  represents  the  left -hand-side  of  the  equilibrium  Eq.  (4). 

The  left-hand -side  of  the  compatibility  equation  is  obtained  from  Eq.  (9)  by 

the  replacement  of  Ag.  by  Ag.  and  A4-  by  A^ . .  This  substitution  is  possible 

because  e2  in  the  compatibility  equation  plays  the  same  role  as  in  the  equilibrium 

equation  and  e  plays  the  same  role  as  M„  .  The  terms  not  involving  A. .  ,  denoted 
1  *  l.l 

f  in  Eq,  (4)  are  given  in  Eqs.  (5).  In  the  development  of  Eq.  (9)  terms  of  higher  order 

than  quadratic  in  the  dependent  variables  have  been  dropped.  Also,  terms  involving 

2 

the  product  of  a  load  and  0  have  been  dropped. 

Boundary  Conditions  ami  Compatibility  Conditions 

Figure  3  shows  a  portion  of  a  shell  meridian  with  a  discontinuity  B-A.  There 
may  be  a  ring  associated  with  this  discontinuity.  The  shear  center  of  the  ring  is 
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located  at  point  C.  External  forces  and  moment  V,  H,  and  M  arc  applied  to  the 

ring  axis  of  shear  centers.  The  shell  ends  A  and  B  also  exert  forces  and  moments 

V  ,  H  ,  and  M  on  the  ring,  those  at  A  denoted  by  a  superscript  plus  and  those 
s  s  s 

at  B  by  a  superscript  minus.  There  are  five  compatibility  conditions  which  must  be 

satisfied  at  the  discontinuity:  compatibility  of  the  two  force  components  and  the 

moment,  and  compatibility  of  the  meridional  rotation  /3  and  horizontal  displacement 

u„  • 

Hr 

If  the  ring  centroid  coincides  with  the  shear  center,  the  horizontal  force 
compatibility  equation  is 


BlluHr  +  B12V  =  -  r/Ts  +  *+/rs  +  H 


(10) 


and  the  moment  compatibility  equation  is 


B12  uIIr  +  B22  V  = 


-  r‘Ms/rs  +  r+Ms/rs  +  M/rs  -  (rVs)  +  flej/rj 


+  (rVB)+(ei  +  ^e2^rs  +  <Me2  ~  ^ei)/r 


+  4>+  (c+  -  ft  e  J  )// 


(11) 


Vertical  force  compatibility  is  represented  simply  by 

<rVg)+  -  (rVsf  +  rgV  (12) 

Compatibility  of  meridional  rotation  and  horizontal  displacement  are  represented  by 
the  following  equations: 


Hr 


"  UH  + 


=  4  +  *+4 


P+  =  fi' 


(13) 

(14) 
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The  quantities  B^  ,  B12  ,  and  B22  in  Eqs.  (10)  and  (11)  are  given  by 
B11  =  EA/rs  +  EIy/ra  "  H^s 


B 


12 


-  El  /r 
xy  s 


22 


EVrs  +  ~T 
Ar 

s 


7ra  -  ^+/r  -  11} 


(15) 


where  E,  A,  I  ,  and  I  are  the  modulus,  cross-section  area,  and  moments  of 
inertia  of  the  ring  about  axes  through  the  shear  center  (see  Fig.  3.)  It  is  clear  that 
Eqs.  (10)  through  (15)  can  be  specialized  for  the  cases  where  the  "discontinuity”  (ring) 
occurs  at  either  end  of  the  shell.  The  displacement  continuity  conditions  Eqs.  (13) 
and  (14)  no  longer  apply  in  these  cases.and  the  appropriate  shell  forces  {either  ( )+  or 
( )  ]  are  set  equal  to  zero,  depending  on  which  end  of  the  shell  is  being  considered. 

The  statement  of  the  nonlinear  boundary  value  problem  is  now  complete.  The 
governing  equations  have  the  form  of  Eq.  (4),  and  the  boundary  conditions  (and 
compatibility  conditions)  are  given  by  Eqs.  (10)  through  (14). 

Solution  of  the  Equations 

All  of  the  equations,  including  the  boundary  conditions,  can  be  written  in  the 

form: 


C +  C2V  +  C3>p  +  C4/3"  +  C5P  +  Cg/3  +  C'jVP  +  Cg^ 

+  C9«i’  +  .  C,/  *  C12  =  0 

(16) 
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The  method  of  finite  differences  is  used  in  the  solution  of  the  equations.  Variable 
station  spacing  is  allowed  for  in  the  analysis,  but  is  not  included  in  the  computer 
program.  Three- point  finite  difference  formulas  are  used: 


*i  =  Vi+1  +  Vi  +  Vi-1 


=  Vi+1 +  Vi  +  Vi-1 


where 


(17) 


C1  =  hi/lhi+l(hi  +  hi+1)J 

c2  =  1/h.  -  !/hi+1 

C3  =  -  hi+l/fhi  <hi  +  hi+1)J 


C4  "  2/<hi+l<hi  +  hi-i  i> I 

C5  =  "  2/\hi+i  U»> 

c6  =  2/(hi(hi  + 


and  hi  ,  hi  +  1  are  the  meridional  intervals  on  either  side  of  the  i1*1  mesh  (joint.  The 
quantities  /?'  and  p"  are  given  by  formulas  analogous  to  Eq.  (17).  Eq.  (16)  in 
finite-difference  form  is 


3  3 


AU*1-1  *  A2lVl  +  *31*1  +  Vl  +  A5l<i*l  +  AeA«  +  2  2  Blk  Wl-*.k 


♦  2  BjVl-2+J  '  Fi 


j  =  l 


jk  i-2+j 

j  =  l  k  =  1 

i  =  1,  2,  ...  K  (10) 
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The  index  i  denotes  mesh  point  number  and  K  denotes  the  total  number  of  mesh 
points  in  the  domain.  Thore  are  two  unknowns  and  PJ  for  each  mesh  point  nnd 
two  equations  (equilibrium  and  compatibility)  for  each  mesh  point  (see  Fig.  4(a)) .  In 
addition,  fictitious  mesh  points  are  introduced  at  either  end  of  the  meridian  and  at 
meridional  discontinuities  in  order  that  first  derivatives  of  ip  and  /?  can  be  specified 
at  the  end  points  of  the  segment  [see  Fig.  4(b)) .  Equations  corresponding  to  these 
points  are  the  boundary  conditions  or  matching  conditions.  Hence,  there  is  a  set  of 
2K  +  4(d  +  l)  nonlinear  simultaneous  algebraic  equations,  where  d  is  the  number  of 
meridional  discontinuities. 

These  equations  are  solved  by  use  of  the  Newton -Raphson  method.  Thurston 
(Ref.  11)  has  discussed  this  method  as  applied  to  various  nonlinear  problems  involving 
many  degrees  of  freedom .  As  applied  in  this  case,  the  Newton-Raphson  procedure  is 
as  follows:  Eqs.  (19)  are  linearized  through  substitution  of 


,n  +  l  ,n  .  .  , 
=  ^  + 


0n  +  1  =  0n  +  A/3. 
Ki  *1  l 


(20) 


where  n  refers  to  "iteration  number"  and  and  A/3,  are  correction  addends  for 
the  current  iteration.  Equation  (19)  is  written  with  A^.  and  A/3,  considered  as  the 
unknowns  and  with  terms  involving  only  >1?  and  /3*1  removed  to  the  right-hand -side. 
Quadratic  terms  in  Ai//.  and  A0j  are  dropped .  The  typical  equation  for  the  current 
iteration  now  appears  as  follows : 

AUA*l-l  +  A2iA*i-l  +  A*3iA*i+  AJl^i  +  A5iA*i+l  +  A6i  A^i+1 

-  Fi  -  All<W  -  A2iCl  -  A3l*?  -  A4L^  -  A5i*M  -  AgAi 


~  2  2  Bjk  ^i-2+J  ^i-2+k  S  Bj^i  ^i-2+j 
j  =  l  1  =  1  j  =  l 

12 
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The  A^.  ,  m  =  1,  6  contain  4k  and  p*}  from  the  previous  iteration.  The  terms  on 
the  right-hand -side  of  Eq.  (21)  are  all  known.  The  simultaneous  linear  set  is  solved 
for  the  vector  of  correction  addends  A ip.  and  A/3.  ,  and  these  are  added  to  the  t'n 
and  /3.  to  give  4? +  and  /3^ +  ^ .  Iterations  continue  until  satisfactory  convergence 
has  been  obtained . 

In  the  Newton-Raphson  method  a  starting  vector  4°  ,  /?°  must  be  assigned.  If 
this  vector  is  zero,  the  solution  for  the  correction  addends,  and  iience  for  4'}  ,  P*  is 

i  i 

the  solution  of  the  linearized  equations.  Problems  that  are  essentially  nonlinear, 
such  as  that  of  a  spherical  shell  submitted  to  a  point  load  at  the  apex,  arc  solved  in 
the  following  way:  Calculations  begin  for  a  small  value  of  the  load,  for  which  the 
linear  theory  is  reasonably  accurate.  Previous  solutions  4k  ,  /?”  are  then  used  as 
starting  vectors  for  each  new  value  of  the  load. 

In  order  that  the  linear  system  represented  by  Eq.  (21)  be  solved  with  a 
reasonable  amount  of  computer  time,  one  must  arrange  the  equations  so  that  the 
matrix  of  coefficients  is  strongly  banded  about  the  main  diagonal.  If  alternate  equa¬ 
tions  are  equilibrium  and  compatibility  equations  and  if  alternate  dependent  variables 
are  A^.  and  A/3.  ,  the  matrix  shown  schematically  in  Fig.  4(a)  will  result.  This 
matrix  corresponds  to  a  shell  with  one  meridional  discontinuity.  The  meridian  with 
mesh  points  is  shown  in  Fig.  4(b).  Each  mesh  point  "generates"  two  columns  and  two 
rows  in  the  coefficient  matrix.  The  matching  conditions  labeled  "Horizontal  Force 
Compatibility,"  etc.  are  given  by  Eqs.  (10),  (11),  (13)  and  (14),  respectively.  All 
elements  not  enclosed  in  a  box  are  zero. 
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Numerous  test  cases  have  been  run  to  check  the  portion  of  the  computer  program 
BOSOR  2  which  calculates  the  axisymmetric  state  of  stress.  The  number  of  mesh 
points  was  varied  to  ensure  convergence.  For  the  test  cases  the  stresses  converge  to 
within  a  few  percent  of  the  exact  values  with  less  than  100  mesh  points.  Two  to  four 
iterations  are  usually  required  for  convergence  to  the  solution  of  the  governing  non¬ 
linear  equations.  Figures  5a  and  5b  show  the  results  of  one  such  test  case.  A  shell 
of  the  geometry  shown  in  Fig.  5a  was  analyzed  by  Stricklin,  et  al  (Ref.  12),  who  used 
the  matrix  displacement  method  in  a  nonlinear  analysis.  In  this  application  of 
BOSOR  2,  95  mesh  points  were  used .  Three  iterations  of  the  Newton-Raphson 
procedure  were  required  before  successive  solutions  were  within  0.1%  of  each  other. 
The  total  computer  time  for  execution  was  3.15  seconds  on  the  Univac  1108. 

Figures  6a  and  6b  show  the  meridional  bending  moments  generated  in  piessurizcd 
steel  hemisphere  with  a  mismatch.  Moments  resulting  from  internal  pressures  of 
100  psi  and  10, 000  psi  are  given  in  Fig.  6a.  The  effect  of  internal  pressure  on  the 
moment  distribution  is  to  shorten  the  "decay  length"  of  the  discontinuity  stress.  This 
is  in  agreement  with  results  obtained  by  other  investigators.  Figure  6b  shows  that 
external  pressure  has  the  opposite  effect:  the  discontinuity  stress  "boundary  layer" 
is  increased  in  length  due  to  the  pressure .  While  the  nonlinear  effect  of  pressure  on 
the  distribution  of  discontinuity  stresses  is  well  known  to  many  specialists  in  shell 
theory,  it  is  rarely  if  ever  accounted  for  in  the  design  of  shell  structures.  This 
effect  could  be  important,  for  example,  in  the  selection  of  positions  of  welds  in 
composite  shell  structures. 
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During  the  manufacture  of  shell  structures  certain  unplanned  variations  in 
geometry  are  likely  to  occur.  Figure  7a  shows  an  example  in  which  the  geometry 
near  the  viewport  of  a  pressure  hull  for  a  deep  submex'sible  vehicle  changed  during 
assembly  due  to  the  weld  on  the  meridian.  The  portion  of  the  hull  in  the  neighborhood 
of  the  apex  moved  inward  slightly  as  shown  by  the  dotted  lines.  The  measured  sinkngc 
is  shown  in  Fig.  7b.  A  curve -fitting  program  (Ref.  13)  was  applied  to  the  measure¬ 
ments  to  calculate  the  variation  in  meridional  curvature  which  is  shown  in  Fig.  7c. 

Two  curves  are  given  in  Fig.  7c,  one  which  corresponds  to  curvatures  calculated  with 
12  input  points  chosen  from  Fig.  7b  and  one  which  corresponds  to  curvatures  calculated 
with  36  input  points  chosen  from  Fig.  7b.  The  nominal  curvature  is  also  shown  in 
Fig.  7c.  The  structure  shown  in  Fig.  7a  was  analyzed  on  the  computer  as  a  shell  of 
two  segments:  the  segment  nearest  the  apex  being  an  "out-of-round"  spherical  shell 
of  variable  thickness,  and  the  other  segment  being  a  perfect  spherical  segment  of 
nominal  curvature  and  constant  thickness  extending  to  a  -  90  deg.  The  thickness 
distribution  is  taken  from  Ref.  14,  and  the  edge  of  the  shell  at  the  viewport  is  assumed 
to  be  free.  The  meridional  bending  moments  for  three  cases  are  shown  in  Fig.  7d. 

In  all  cases  the  shell  is  submitted  to  uniform  external  pressure  of  1.0  psi.  The  bonding 
moment  distribution  for  the  perfect  shell  agrees  with  that  calculated  in  Ref.  14.  This 
moment  distribution  is  compared  with  moments  for  the  imperfect  shell  in  which  12 
and  36  input  points  are  used  for  calculation  of  the  geometry  of  the  variable  thickness 
segment.  Two  interesting  observations  might  be  made  with  regard  to  the  results : 
the  stress  distribution  is  rather  drastically  affected  by  the  weld  sinkage  (the  maximum 
stress,  which  occurs  near  the  weld  at  the  outer  fiber  being  increased  by  almost  a 
factor  of  two;  and  the  moment  distributions  corresponding  to  the  two  different 
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geometries  (Fig.  7c)  calculated  by  the  curve-fitting  routine  are  very  similar,  even 
though  the  input  meridional  curvatures  are  very  dissimilar.  Apparently,  the  shell 
behavior  is  related  to  a  "smoothed"  or  averaged  meridional  curvature  distribution, 
rather  than  to  the  rapidly  changing  variation  of  curvature  shown  in  Fig.  7c  as  a  dotted 
line.  The  probable  explanation  is  that  the  arc  length  over  which  the  rapid  curvature 
changes  occur  is  small  compared  to  the  attenuation  length  or  "boundary  layer"  length 
of  the  discontinuity  stresses  for  this  particular  shell. 

Figure  8  gives  the  hoop  stress  resultants  for  pressurized  shallow  spherical 
shells  with  edge  rings  located  as  shown.  Edge  moments  are  also  applied  as  shown. 

The  vertical  reaction  coincides  with  the  ring  centroid,  which  is  located  distances 
and  from  the  edge  of  the  shell.  The  pressures  are  0.619  psi  and  0.357  psi  for 
the  two  cases.  It  is  seen  that  small  differences  in  the  geometry  of  the  structure 
cause  large  differences  in  the  stress  distributions.  The  ring  eccentricities  must  be 
accounted  for  in  the  analysis  of  such  shell  structures.  As  an  aid  to  the  designer  the 
computer  program  can  be  used  to  find  the  most  favorable  location  of  the  edge  ring. 

Figure  9  shows  the  maximum  meridional  stress  in  the  outer  fiber  of  an  internally 
pressurized  aluminum  cryogenic  tank.  The  geometry  is  shown  at  the  top  and  bottom 
of  the  figure.  The  inner  surface  from  the  centerline  to  the  point  D  is  an  ellipsoid 
and  from  D  to  the  plane  of  symmetry  is  a  cylinder.  The  thickness  varies  linearly 
between  the  stations  where  it  is  called  out.  The  large  increase  in  thickness  near 
the  juncture  of  the  cylinder  and  the  ellipsoid  was  included  in  the  design  in  order  to 
compensate  for  the  weakening  effect  of  a  weld  there.  The  outer  fiber  meridional 
stress  is  maximum  at  point  C.  As  seen  from  Fig.  9  linear  theory  yields  93  ksi  for 
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this  maximum,  and  nonlinear  theory  yields  74  ksi.  The  results  for  linear  theory  arc 
in  good  agreement  with  a  solution  obtained  through  use  of  a  computer  program  written 
by  Kalnins  (Ref.  6). 

Figure  10  shows  SC  4020  computer  plots  of  stress  resultants  for  the  cryogenic 
tank.  The  plot  routine  was  written  by  Dyche  (Ref.  15).  Various  configurations  were 
analyzed.  Figures  10b  and  lOe  correspond  to  a  tank  in  which  the  reference  surface  is 
an  ellipsoid  joined  to  a  cylinder  with  the  thickness  distributed  equally  on  both  sides  of 
the  reference  surface.  Figures  10a,  c,  d,  f,  and  g  correspond  to  a  tank  in  which  the 
same  reference  surface  geometry  is  considered  to  be  the  inner  surface  of  the  shell. 

As  seen  from  a  comparison  of  Figs.  lOe  and  lOf,  the  large  bending  stress  at  point  C 
in  Fig.  9  disappears  if  the  shell  is  manufactured  such  that  the  reference  surface  is 
the  middle  surface.. 

Figures  lOd  and  lOg  correspond  to  1/1000  of  the  design  pressure  p  =  70.3  psi, 
and  thus  approximate  the  linear  solution.  Note  that  the  hoop  stress  resultant  calculated 
from  linear  theory  is  negative  in  a  small  region,  giving  an  erroneous  indication  of  the 
possibility  of  buckling  at  some  higher  pressure.  Also,  the  hoop  stress  resultant  at  the 
plane  of  symmetry  would  be  over-estimated  by  about  25%  if  linear  theory  were  used 
for  calculation  of  stresses  at  the  design  pressure.  The  bending  stresses  predicted 
by  the  linear  theory  are  about  twice  those  predicted  by  the  nonlinear  theory. 

Figure  11a  shows  a  portion  of  a  spherical  shell  with  a  local  area  in  which  the 
radius  of  curvature  is  greater  than  the  nominal  radius.  This  "flat"  area  represents 
an  imperfection  in  the  shell,  which  is  submitted  to  uniform  external  pressure. 

Figures  lib  and  11c  give  the  state  of  the  shell  just  before  it  collapses  in  an 
axisymmetric  mode.  Collapse  loads  are  calculated  as  described  in  Ref.  2. 


17 


LOCKHEED  PALO  ALTO  RESEARCH  LABORATORY 

tOCkMttb  11111111  i  IMCt  COMPANY 
»  oto«r  •ivitioN  o»  loimiit  iuckm  co«*oa*f>nN 


In  the  plots  the  origin  corresponds  to  the  apex  of  the  shell.  The  shell  is  treated  as  if 
consisting  of  three  segments:  the  "flat"  spot,  a  region  approximately  one  "boundary 
layer"  in  length  adjacent  to  the  flat  spot,  and  the  remainder  of  the  shell.  The 
meridional  station  which  corresponds  to  the  edge  of  the  "flat?*  spot  is  indicated  by 
heavy  vertical  lines  in  Figs,  lib  and  11c.  The  dotted  lines  in  the  plots  of  UH  and  UV 
represent  the  membrane  solution  extrapolated  to  the  apex  of  the  shell.  The  plots 
correspond  to  a  load  within  0.1%  of  the  collapse  load,  which  is  28. 7%  of  the  classical 
buckling  load  for  a  complete  sphere  of  the  nominal  radius.  In  Ref.  16,  a  collapse 
load  of  28.0%  of  the  classical  buckling  load  was  found  in  an  analysis  in  which  the 
entire  shell  was  treated  as  one  segment.  It  is  seen  that  the  maximum  inward  dis¬ 
placement  (UV  at  the  apex  =  0. 85)  is  considerably  greater  than  the  maximum  value 
of  0.7  shown  in  Fig.  2  of  Ref.  16.  Near  the  collapse  load  the  displacement  increases 
rapidly  with  small  increments  in  the  pressure.  This  case  is  included  to  illustrate 
the  advantage  of  analyzing  a  single  shell  structure  in  segments.  The  input  data  is 
simpler,  and  fewer  mesh  points  are  required  for  convergence  to  within  a  prescribed 
error. 
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Fig.  1  Typical  "Practical"  Shell  Structure 
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Fig.  4(a)  Coefficient  Matrix 
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Fig.  4(b)  Meridian  with  Mesh  Points 
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CIRCUMFERENTIAL  STRESS 


Fig. 


<b) 


5  Meridional  and  Circumferential  Stresses  in  a  Pressure  Vessel 
With  Internal  Pressure.  Comparison  with  the  results  of  Ref.  12 
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Fig.  6  Effect  of  (a)  Internal  and  (b)  External  Pressure  on  the  Discontinuity  Bending  Moment 
Sphere  With  a  Mismatch 


(a)  GEOMETRY 


WELD 


Fig.  7  Effect  of  Weld  Sinkage  on  the  Bending  Moment  Near  a  Viewport 
of  a  Pressure  Hull 


29 


LOCKHEED  PALO  ALTO  RESEARCH  LABORATORY 

IOCIMID  Mlltllll  A  I  H  C  I  COMPANY 
k  oio«»  ftivttto*  iortHif.. 


HOOP  COMPRESSION,  N?n  -  LBS/IN 


WANG  SPHERICAL  CAP, MO  =0.8 


0  2  4  6  8  10 

ARC  LENGTH,  s 


Fig.  8  Hoop  Compression  for  Cases  1  and  2 
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Fig.  9  Cryogenic  Tank  and  Maximum  Meridional  Stress  From  Linear  and  Nonlinear 
Theories 
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Fig.  10  Stress  Resultants  in  Cryogenic  Tank 
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Fig.  11  Axisymmetric  State  of  Shell  With  "Flat"  Spot  at  Collapse 

(b)  Radial  Displacement,  Axial  Displacement,  and  Meridional 
Stress  Resultant 
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